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SUMMARY 

An exact treatment of the problem of finding the incompressible, 
inviscid two-dimensional flow around two cascades in tandem is pre- 
sented. The analysis includes solutions of both the direct and the 
inverse problems. These problems are solved by conformally mapping 
the tandem cascade onto the region between two concentric circles in 
which region there are suitably placed flow singularities. Formulas 
for the velocity and the potential in the annular region are presented 
in a closed form by rneems of elliptic functions. The eqixations are 
presented in a form suitable for computation. 


INTRODUCTION 

A theory of the two-dimensional flow of an incompressible inviscid 
fluid through two cascades in tandem, which was developed at the NACA 
Lewis laboratory, is presented herein. Methods are developed for the 
solution of both the direct and the inverse problems for a tandem 
cascade. For a. single cascade, there have been several solutions of 
these problems. A survey of such problems is presented in reference 1. 

The results of this analysis can be applied to a variety of 
problems that occur in compressors. The amount of ttiming that pan 
be acconiplished by a single row of blades without separation of the 
blade boundary layer appears to be limited. Criterions for the amount 
of circulation attainable around a blade of a cascade without separa- 
tion of the boundary layer, with the attendcuit increased losses, 
are given in reference 2. If an amount of turning greater than that 
which can be obtained with a single blade row without boundary-layer 
separation is desired, this turning may be accomplished by the use of 
two rows of blades, which are so designed that blade boundary -layer 
separation does not tsike place. 

Another situation to which the results of this analysis are 
applicable is the following: In some compressors, there are additional 

rows of turning vanes beyond the last row of stator blades. These rows 
of tTirning vanes form a configuration to which the present analysis can 
be applied as a two-dimensional approximation. In particular, indica- 
tions of flow for angles of attack other than the design angle of 
attack can’be obtained. 
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A tandem cascade is composed of two cascades arranged approximately 
as shown in figure l(a) . The method -employed herein for the study of 
this problem is the method of conformal transformation of the tandem 
cascade into a standard doubly connected region for which the flow is 
known. Various stages of the transformation are sketched in figure 1 . 

In many respects, the problem of the tandem cascade is similar to that 
of the biplane. As is apparent in figinre l(b), the first stage of the 
transformation renders the tandem cascade into what is essentially a 
two-dimensional biplane. The standard doubly connected region chosen 
is the same as that used in reference 3 for the biplane problem; that 
is, the region contained between two concentric circles. 


GEHEEAL EQUATIONS FOR FLOW THROUGH TANDEM CASCADE 

In order to determine certain relations between various flow para- 
meters of a tandem cascade as shown in figure l(a), the following 
procedure is used. The lines P1P4 and P2P3 are located far down- 
stream and far upstream of the tandem cascade, respectively, (fig. 1(a)). 
The streamlines P^Pg and P3P4 are located exactly one cascade 
spacing apart; P1P4 and P2P3 are parallel to the cascade axes. 

From the continuity condition, the net flow from the contour 
P1P2P3P4P1 1 (a)) niast be zero. Because P1P2 and P3P4 are 

streamlines, the flow across them is zero. The flow through P2P3 
must therefore equal the flow through P1P4. From this relation 

tVi cos Xj- = tV2 cos X2 

where t is the cascade spacing (the same for both cascades) or 


Vi cos X^ = Vg cos X2 (1) 

where and Xq are the limiting values of magnitude and direction 

of the velocity far upstream of the tandem cascade, respectively; and 
Vg and A2 limiting values of the same quantities downstream 

of the tandem cascade. (The main symbols used herein are defined in 
appendix A. ) 


The circulation around the blades enclosed within the contour is 
given by 


r 


Pl^’2P3Vl 


u ds 


( 2 ) 
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The integrals €Q.ong P 1 P 2 and P 3 P 4 > in equation (2), are equal in 
magnitude and opposite in sign. This integral therefore reduces to 


or 


/■ 


P1P2P3P4P1 “ 


ds 


t sin ~ ^2 sin ^2) 


r 

t 


sin X 2 . ■ ^2 sin }^2 


( 3 ) 


In figure (2), Vm and refer, respectively, to the magnitude and 

the direction of the vector mean of the upstream and downstream velocities. 
By using complex-number notation, the mean velocity can he written as 


iXn iXp 


( 4 ) 


Also hy using the continuity condition (equation (l)), equation (3) 
can he written: 


iT 

2t 




2 


Solution of equations (4) and (5) together gives 


( 5 ) 


r 

Vn e = V e +i — 
''1'=' +^2t 


( 6 ) 


Voe^^ = V e^^^-i— 
''2® '^m® -^2t 


( 7 ) 
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FLOW IN BEGION BOUNDED BY CONCENTRIC CIRCLES 

Derivation of equations for velocity and potential . - The veloci- 
ties on the ho\mdaries of an annvilar region, in which there are placed 
two complex sources, must he fo\uid. By means of successive reflections 
of the singularities in the two boundaries, ein expression will be 
obtained for these velocities. The details of the process of reflec- 
tion will be presented for only one of the singularities inasmuch as 
the process will be the same for both. In figure 3, A and C are, 
respectively, the radii of the iruaer and the outer circular boundaries 
of the annular region, which is so located that the centers of these 
circles are at the origin. The original singularity is located at B. 

For convenience, the real (source) and the imaginary (vortex) 
parts of the singularity are considered separately. First to be 
considered is the reflection of the source located at B. The bound- 
ary condition on the surface of radius A is that there should be no 
conponent of velocity normal to this surface. In order that this 
condition be fulfilled, a like source must be added at a point inverse 
to B“ with respect to the circle A and a source of opposite sign 
must be added at the center of the circle. The boundary condition on 
surface C is the same as on surface A. It is apparent that the 
first reflected sovirce now violates the boundary condition on the 
outer boundary. This image singularity must therefore be reflected 
in the outer boundary, which again necessitates a reflection in 
the inner boxmdary and so on indefinitely. If the reflection 
process is considered to begin with a reflection in the outer boundary, 
auiother series of images in different positions results. An additional 
singularity has to be placed at the center of the circle for each step 
of the reflection process. The total flow across boundaries A and C 
must be zero. The real parts of the two original singularities at 
B and B’ are therefore equal in magnitude and opposite in sign. The 
reflection is carried out in the same manner for the other original 
singul€irity located at B*. 

The reflection scheme for the vortex is similar to the one for the 
source except that on each reflection the sign of the reflected vortex 
is changed and the singularity added at the center of the circle has a 
sign opposite to that of the reflection. The magnitudes of the distance 
from the origin to the reflections are shown in figure 3. 

Let 

X source strength of singularity at B 

-Y vortex strength of singularity at B 
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X' source strength of singularity at B’ 

-Y’ vortex strength of singularity at B' 

As previously discussed, the following relation must exist: 

X = -X' (8) 

The potential at a point z due to a source of strength X 
located at a point zq is^ if the potential is to he zero, at a 
point zj_, 




z-zq 

Zi-Zq 


For the vortex -Y located at zq, the potential at z is, if it is 
desired that the potential he zero at a point z^. 


2jt 


loge 


z-zq 

zi-zq 


The potentieuL due to the sovurce at B and its reflections can he 
written (see fig. 3 for locations of reflections) 


Wx(z) 



r- 00 

\ \ 

1 A^ 

X 


A B 

2rt 


z 1 A2n 


Z- / 

1 

L ^ B c2(a-l)J 


+ 



1 


z 


A 



1 ~ 
B A^(q-I) 

1 

B 



( 9 ) 
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The point A vas selected as the point for which the potential is zero. 
The potential due to the vortex at B and its reflections is 





( 10 ) 


There are similar expressions for the ccmplex potentisil that result 
from the source and the vortex at B' . 

For the total potential due to hoth singularities 


W(z) 
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These series are absolutely convergent in a closed region that excludes 
the singularities B and B’. Differentiation of equation (ll) yields 
the conjugate velocity 



Velocity and potential in closed form . - In order to evaluate equa- 
tions ^11) and (12) more readily, the annular region in the z -plane is 
transfotrmed into a rectangle in the ¥-plane. The transformation used is 


z 
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By introducing the notation 

V(w) = V[z(w)] = V(e‘'^) 

equation (U) can he written as a fianction of v 



(13) 


The potential can also he written in closed form as a fiinction of w. 
W(w) = W(eW) = 2+iI logg 0 ^ (w - h) + iK' + 

logeS ^ (w+h+2a)+iK^ + logeS ^(w-h')+iK> + 

loggS (w+h'+2a) + iK' + Aiw (14) 

where Aj_ is an imaginary constant. Equations (13) and (14) are 
derived in appendix B. The symbols for elliptic functions are those 
defined in reference 4. 
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SOLUTION OF INVERSE PROBLEM 

The inverse problem, as considered herein, is to find a pair of 
profile shapes that when suitably arranged in a tandem cascade will 
have a prescribed turning and surface velocity distribution. The 
region external to the tandem cascade is to be conformally transformed 
into a standard region through which the flow is known. The annular 
region between two concentric circles with suitably placed flow sin- 
gularities is chosen as the steihdard region: , (l) The general form of 
the mapping function together with certain conditions upon it will be 
found. (2) A method of computing the values of the mapping function 
along the boundaries of the annular region is given. (3) Because the 
assigned conditions may not give a closed profile, a method is given 
by means of which the initially assigned velocity distributions may 
be modified to give closed profiles. 

Derivation of mapping function . - The a-plane is the plane of 
the original tandem cascade. In order to transform the tandem cascade 
into the standard region, the following series of transformations is 
used: 

( 1 ) loge i = ^ (15) 

Because of the period properties of the exponential function, equa- 
tion (15) transforms the tandem cascade into two closed figures in the 
|-plane (fig. 1(b)). Under this transformation, points for which 
Re(o) is very large positively are transformed into points that are 
a large distance from the origin in the | -plane. Points having a very 
large real part that is negative are transformed into points near the 
origin. Therefore, 

lim I (a)= 0 B 

Re(o)-* + oo (16) 

lim I (a)=0 

Re(a)-»-po (17) 

(2) According to reference 5, it is possible to transform conform- 
ally the region exterior to the closed figures in the | -plane into the 
area external to two circles. It is convenient to write this trans- 
formation in the following form: 
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One circle, the correspondence of one point thereon, and in addition, 
the location of the center of the second circle cein he specified. If 
the point | = ® is transformed into some finite point of the 
tl -plane, a bilinear tremsformation can be used to transform this 
finite point into the point at infinity in the t 2 -plane. Under a 
bilinear transformation, circles are preserved as circles. A trans- 
lation and a rotation can then be used to place the circles as shown 
in figure l(c) where the point at infinity has been transf onned into 
the point at infinity in this plane, the t -plane, and the centers of 
the circles lie on the axis of real niambers, equidistant from the axis 
of imaginary numbers. This entire series of transformations can be 
combined arid written as 




(18) 


( 3 ) Circles C and D of the J-plane can be transformed into con- 
centric circles in the z-plane by 



(19) 


where ki and kg are points such that they are inverse points with 
respect to either circle. Therefore, 

(ki-k) (kg-k) = 

(ki+k) (k2+k) = 

where Rq and % eire the radii of the two circles in the J -Plane. 
These two equations may be solved for ki and kg. 

If equation (19) is solved for t ia terms of z, 

kgz+ki 
^ ■ z-1 


dt _ ~(k2+ki) 
dz (z_i )2 


From which 
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and therefore, 


( 20 ) 


^ ^ d| ^ ^ .e^(^)(k2 + ki) 
dz dt dz (2 _ 1)2 

When the notation f(z) is used for f [J(zj], equation (20) becomes 

d| _ ) 

(z - 1)2 

If the point to which | = 0 is transformed is denoted hy Re^^ 

I = (23) 

where f(z) is an entire function and 


^ " ' ( ' ^:t)2 - 1) F'(z) + (Rei^ - 1)] 


d| 


( 21 ) 


(23) 


Equating these two expressions for and solving for F’(z) gives 

dz 


F'(z) = 


-(k2 + ki) e^(^) - F(z) - (Reir. i) 


(z - Re^^)(z - l) 
This equation can he rewritten as 


_ 1^3 (1 - e8(^)) 


F«(z) = 


(z - Re^^) (z - 1) 


(24) 


k-r = 1 - Re^^ 


where 
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and 


>g(z) =. (ki + k 2 ) e^(^) ~ 
■ 1 - Re^r 


(25) 


Because F(z) and therefore F’(z) are regular in the anniolar region, 
1 _ eS(z) must have zeros at 1 and Re^lf*. Therefore 


g(l) = 0 
g(Rei‘>^) = 0 


Because 


therefore. 


loge I = 


2jccy 


(26) 

(27) 

(15) 


' 2jia 

^ ^ e' 

^ 2jt 

Therefore, hy using equations (22) to (25) and (28) 

dz d| dz 2n(z - l)(z - Re^l^) 


(28) 


(29) 


The conditions on g(z) may be more easily applied if the function 
Q(z) is introduced such that 


g(z) = Q(z) h(z) 


(30) 


vhere h(z) is analytic in the annular region and Q(1) and Q(Re^^) 
are both zero. In addition, Q(z) must be read on both boundaries. An 
example of a function siiitable for Q is presented in appendix C. 


Method of cong)utation of mapping function . - The cr-plane and the 
z -plane are related by a conformal transformation 


U((t) = W(z) 
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where 

U(o) potential function in tandem-cascade plane 

W(z) potential function in annular-region plane 

Differentiation of the preceding equation yields; 

dU(g) ^ ^ dW(z) 
do dz “ dz 


or 

u(a) 

Let u(z) = u(cr(z) ) . 

Equation (31) can he rewritten as 

dz u(z) 


From equation (29) 


(31) 


(31a) 


^ = ± ee(z) 

dz 2rt 


^ 1 z - Re^^ 


(32) 


From eqviations (16) and (17) and the discussion preceding equation (l8), 

iim . . z=l ^ . 

Re(a)-»’+<» 

Iim zpReiT 
Re(a)-*.-oo 


From equations (6), (7), and (31a), 


lA, 


m 


ir 


Iim u(z)=Vnie + ^ 
z-*-l 2t 

^m 

lim u(z)=Vine - ^ 
z-Rgir 
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so thalb the residues of V(z) at 1 and Re^"** are, respectively, 

t - |T) 

t .. g) 

If the points 1 and Re^^ are. identified with the points B and B* 
of equations (ll) and (l2), the source and vortex strengths are 

X = -tVm cos Xjjj 

Y = ^ + tVm sin 

X' = tVm cos Xm 

Y‘ = - tVm sin Xm 

Next, an annular region such that the ranges of potential are cor- 
rectly matched to the assigned ranges of potential on the two blade 
shapes must be foimd. The following scheme may be used; Let Re^Y 
equal -1. The radii of the inner and outer boundaries are left to be 
determined in a manner such that the ranges of potential have the proper 
values . 

A certain amount of control over the relative positions of the two 
cascades can be achieved by specifying the difference in the stream 
function for the two blades and the potential difference between the 
tail stagnation point on one blade and the tail stagnation point on the 
other blade. These additional specifications increase the difficulty 
of determining the boundary radii and the location of the point Re^Y. 

Once the radii of the boundaries and the locations of the singulari- 
ties have been determined, the correspondence between points on the two 
blades and points on the circular boundary may be found by matching 
points having the same potential (as in reference 6) . From this cor- 
respondence is known' the assigned blade velocity as a function of the 
central angle on the corresponding circle. 
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From equations (29) and (30) 


da -tks eQ(z)^Cz) 
dz ■ 2jt (z - l)(z - ReiT^) 


( 34 ) 


where 


li(z) =2, 

— 00 

= P(x,y) + iq(x,y) 

7(z) ^ ^ 

u(z) ~ dz 

On the boxmdaries of the anntjlar region, Q(z) is real and 


(35) 

(31a) 


Vfz) 


1 


1 


gQ(z)h(z) 


-tk3 

u(z) 


z-1 


z - Re^T* 



2jt 


(36) 


Therefore, on either circuleor houndary 


r 


Re h(z) = ^ 


loge 


s 


z) 


+ loge 


z-1 


+ loge 


Re 


ir 


- loge 


-tkg 

2rt 

(37) 


By use of equation (37), Re (h(z)) can be computed on both circular 
boundaries. In order to compute the fimction conjugate to Re (h(z)), 
use is made of Villat's analogue; for the annular region, of Poisson's 
integral. These equations are derived in reference 7. A form more 
convenient for the present purpose is given in reference 3 (p. 13, 
equations (23) and (24)). 
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if'" 

^A ('P') =^o + § Pc(^) Z (cp - cp«) dCp - 

0 

2jt 

(cp - cp' ) dqj ( 38 ) 




p2jt 

q^(cp') = bo + i J pc(cp) (cp - cp' ) dCp 
0 


1 r 

rt PaM Z (cp - qp' ) dcp 
0 


where 


(39) 


P^(cp) real part of h(z) on inner boundary, as function of central 
angle cp 

Pp(qp) real part of h(z) on outer boundary as function of central 
angle qp 

q^(cp) imaginary part of h(z) on inner boundary as. function of 
central angle cp 

q^(qp) imaginary part of h(z) on outer boundary as function of 
central angle cp 


Zi(u) 


constant^ which may he taken as zero 
H'(u) 




(See reference 4 for definition of the H function.) 


From equation (34) 


-tk, + iq(x,y)l 

dcT = — dz 

2jt (z - l)(z - Re^^) 


( 40 ) 
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Integration of equation (40) gives the required value of a. 

In order to find the position of one cascade relative to the other, 
it is necessary to integrate ^ along a line from one circular 

CjLZ 

■boundary, across the annular region to the other circular "boundary. 

This procedure involves knowing the values of h(z) along this line. 
These values can "be found by means of the following equation, which 
(with different notation) is found in references 6 and 3; 


h(z) = (ao + Ibo) 



p^(cp) Z i log(j) + dcp + 


where 




1 

it 


p^((p) zji iog(^) +cpj dcp 


u 

0 






Pa(«P) dV = i 


2jt 

r> 


Pp(cp) dcp 


0 

0 


0 


(41) 


Modification of assi^ed velocity to obtain closed profiles. - 
The preceding development does not guarantee that the resulting blade 
profiles will be closed curves. In order to modify the Initially 
assigned velocity so that closed blade shapes will result, the follow- 
ing procedure is used; 


The sum of the residues of the derivative of the mapping function 
inside the contour indicated in figure 4(a) is zero; that is, the 
Integral around the contour will be zero. The contour in the a-plane 
corresponding to the contour indicated in figure 4(a) therefore will 
be closed, whereas the portions of the contour corresponding to the 
bounding circles inl^t not be closed cuirves. If this part of the con- 
tour is open, the opening for one blade' will be the negative of the 
opening for the second blade because the integration around the outer 
circle is in the negative direction (fig. 4(a)). 
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The amount hy which the profiles fail to close can he reduced to 
zero hy adding to the derivative of the mapping function, a function 
having a pole at the origin such that the residue of this fxanction at 
its pole is eq^aal to the negative of the amount of opening divided 
hy 2it, 

If the amount of lack of closure is eqiial to m^, the modified 
mapping-function derivative can he written as . 

dg' _ ^ _ V* (z) 

dz “ dz “ 2jtz “ u' (z) 

Because the velocity remains undisturbed in the annular region 

V*(z)=V(z) 


Therefore, the altered velocity is 


^ z V(z) - mo u(z) 


and the altered arc length is given hy 


dg* = 


2jtz 


VCz) 

ufzi 


dz 


(42) 


(43) 


SOLUTION OF DIRECT PROBLEM 

The direct problem is that of finding the flow of fluid past blades 
of known shape. In addition, to the blade profiles, the upstream and 
downstream velocities are also known. The Kutta condition of finite 
velocities at the trailing edge of the blades is used to^ determine the 
blade circvilations . 

As in the section "Solution of Inverse Problem", by the application 
of an exponential transformation, the flow through the tandem cascade 
is transformed into the flow about two closed profiles with a flow 
singularity in the finite part of .the plane; this is essenti al ly the, 
problem of the two-dimensional biplane. A method by which the trans- 
formation from the region external to the two closed shapes to the 
region between two concentric , circles can be con5>uted by means of 
successive approximations to the solutions of a pair of Integral equa- 
tions is given in reference 3. 
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In the biplane problem, the flow singula r ity is a single dipole. 
In the problem tinder discussion, there are two singularities in the 
flow field. In order to determine the location of the addition singu- 
larity, use must be made of a trial -and-err or procedure; an equation 
similar to equation (41) is used to determine the values of the map- 
ping function at points interior to the annular region. From equa- 
tions (l3) and (l4), the velocity and the potential on the annular 
region boundary is Imown. The analysis is completed by transforming 
these velocities back to the original tandem cascade configuration. 


Lewis Fli^t Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 9, 1951. 
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A, C 
a,c 

B, B* 
t>,b* 

m 

f(z),F(z),g(z), 

h(z),Q(z) 

G(w) 

Ii3(w),h4(w) 

niQ 

p(x,y) 

Pa('P) 

Pc(<?>) 

<l(x,y) 

%{<P) 

r 

t 

U(a),W(z) 


APPEaroiX A 
MAIN SYMBOLS 

radii of boundaries in z-plane 
a. = - logg A c 2 logg C 
complex constants 
singularity positions in z-plane 
b s logg B b' s logg B' 
analytic f inaction of ^ 
analytic functions of z 

velocity function in w-plane 
arbitrary entire functions 


residue of function added to achieve closure 
Re [h(z)] where z = x + iy 

Re [h(zj] as function of along jz| = A, z = Ae^*^ 

Re [h(z3 as function of along |z| - C, z = Ce^*^ 

Im [h(z)] where z = x + iy 

Im Cli(z)] as f inaction of <P along |zl - A, z = Ae^*^ 

Im [h(z)] as function of cp along |z| = C, z = Ce^*^ 

singularity position in z-plane (Re^^ = B') 

logg R 

spacing 

potential functions 
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V(z),u(a),u 

velocity fuQctions 

V(w) 

V(w) = V[z(w)J 

u(w) 

u(w) s u[a(w)] 

Vie^^l 

upstream velocity 


downstream velocity 


vector mean velocity 

X,X' 

source strengths 

-Y,-Y’ 

vortex strengths 

r 

circulation 

0-plane 

tandem-cascade plane 

0' -plane 

modified tandem- cascade plane 

z -plane 

annular-region plane 


auxiliary variables used in transformation 


Notation used for elliptic functions is that found in reference 4 
except for Z 2 _(u), which is defined in text following equation (39). 



22 


NACA ra 2393 


APPEiroiX B 

EmiESSION OF VELOCITY AMD POTEMTIAL IM TERMS OF ELLIPTIC FUMCTIONS 

Velocity expressed in closed form . - In order that equation (l2) 
may be more readily evaltiated, the annular region in the z-plane is 
transformed into a rectangle in the v-plane. The transformation used is 

z = e'^ (Bl) 

Let 

V(w) = V[z(v)] = V(eW) 

Then equation (12) can be rewritten as a function of w: 
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Let 

Further let. 

Then, because X 
(fig. 1(e)) 

because 

If 

Then 
and 
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logg A = -a 
loge B = b 
loge C = c 
loge B' = b* 


-jt < Im(w) < rt 
-a < Re(w) < c 

= -X' , for V contained within this rectangle 


V(w+2a+2c) 



(B3) 


gW+2ni ^ gW 

V(w+2ni) as V(w) 


G(w) as V(w) (B4) 

G(w+2ni) = G(w) (B5a) 

G(w+2a+2c) = V(w+2a+2c) 

= e'^ V(w) ^ 2i = G(w) (B5b) 


Therefore G(w) is a doubly periodic or elliptic function and has a 
real period (2a+2c) plus an imaginary period 2ni. 
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Consider one of the terms of G(w) 


cW 


- B. 


p2n 


(B6) 


This ftinction has an expansion in terms of jjw - b + 2n(a + c)] 


-1 + 


= _i + f [v - b + 2n (a + cFI 

I ^ ^ + 2n (a +.C) kl r 

fi2n \ Lk=0 J 


[w - "b + 2n (a + c)]^[^ = 

k=0 


w - b + 2n (a + c) ^ 


-1 


(B7) 


where the are Bernoulli numbers (reference 8, p. 183). The first 

Bernoiilli number Bq is equal to 1. The residue of the fimction of 
equation (B6) is therefore 1. Similarly, the residues of the functions 


=.w 


=.w 


e'^ -B 


e^ - B 


i2n 

j 

\2n 


> and 




pW i 

■ B c2(n-l) 


- 1, and so forth are each 1. 


Consider a parallelogram in the w-plane, 

-irt iS Im(w) < n 
Re b + 2a < Re(w) ^ Re b + 2c 

In_this rectangle G(w) has four poles located at b, - b - 2a, b*, 

- b' - 2a. The sxmi of the residues of G(w) at these poles is 

I 

(X + iY + X - iY + X' + iY' + X' - iY') = 0 (b8) 

because 


X = -X' 
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The principal part of G(v) at each of these poles is, respectively, 

X + iY 
2rt w - h 

A. X - iY 
2rt w + h + 2a 

2^ X' + iY' 

2n w - h' 

1 . X' - iY^ 

2rt w + h? + 2a 

« 

According to reference 4 (p. 474), if 


^r 

Ar,m ” Pr) ^ 

nt=l 

is the principal part of an elliptic f\inction f(z) at its pole Pr> 
then 


n 


f(z) = Ai + ^ 
r=l 


mr n 

m=l ' o-z 


Atz - jtPj. 

«^2 \ 

1 2o>l 

“1/ 


where A^^ is a constant . 
From reference 4 (p. 479), 


(B9) 


H(z) ='>^1 2^3 


-2 




^ = JL 

3 2K 

so that, hy using the plder notation of Jacobi, 


— ' \ — ' , - ^ ni-1 


V V (-1) d® 

f (z) = Ai + ^ (m . 1) ; Ar,m l°Se ^ 


r=l m=l 


K r 

— (z 


- Pr) 


iK’ 

K 


(BIO) 
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For the present case, 


A - X+iY 

Al,l - — 

A X-iY 

2rt 


so that. 


, X'+IY* 

^3,1 - 2n 

. .. X'-IY' 


H,i - 2r~ ■ 

Pi = t 

Pg = _b - 2a 

P3 = V 

p 4 = -h' - 2a 


(x+iY\ 

H' 

— (w - b) 


H' 

— (w + b + 2a) 
coi ^ ' 

\ 2« ) 

H 

— (w - b) 

Oil '■ ' 

^ \ 2n j 

H 

^ (v + b + 2a) 

L“l J 


:'+iY'\ 

2n J 


H' 


H 


Lcoi 


(w - h'^ 


(v . h-) 

\9i j • 


/X’-iY’ 
\ -.211 


H' 


— (w + h' + 2a) 

Oil ' 


(w + h' + 2a) 


(Bll) 


or 




0 ^ 


(w-b)+iK' 




^ (w+h+2a)+iK* + 




^ (w + b' + 2a) + iK'l 


(B12) 
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where 

Z(w) 


When Jacohian elliptic functions are introduced into equation (B12) and 
equation (B4) is used in conjiinction, 





, (13) 

Let 


w = -a 


-ifl H»(w + iK») 

2k H( w + iK' ) 
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(B13) 


Por convenience let 


— Im b + K’ = V 

% 


— Im b' + K' = v' 


— Re b + a = u 


= u' 


— Re b' + a 

o>L 
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Then equation (B13) can he rewritten as 



I 

L L 

iY 

Z(u) 

. if 


' TzCiv) + 

CO 


V sn(iv) sn^(u) cn(iv) dn(iv) 

^ 1 - k2 sn^u sn2(iv) 


1 - k^ sn^(u) sn^ (iv) 


1 - k^ sn^(iv') sn^(u') 


/ ,y _ k^sn(iv») sn(u') cn(u* ) dn(u')~ |l 
1 - k^ sn^(u') sn^(iv') _J^ 


(B14) 


The point w = -a corresponds in the z-plane to the point z = A. At 
this point, the velocity, in order to he tangent to the boundary, can- 

2K 

not have a real component. The coefficient of — in equation (B14) 


OJn 


is an imaginary number. Therefore, must he an imaginary number. 

Re{Aj) = 0 (B15) 


dW(z) _ - 


dz 


= V(z) 


dWM2iL ^ = v(z) 
dw dz ' 


(B16) 

(B17) 


From equation (B1) 


dz „w 
= 


(B18) 
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If in equation (U), z is replaced by- e^, eqiiation (ll) is equiv- 
alent to 


X+iY X+iY 



Consider the first bracketed term of equation (B19) . This 
product has zerps of the order of 1 at points congruent to 
b modulo 2(j(a + c) + krti), where j and k are integers. The 
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fmction (w - 15 ) + iicij has zeros at the same points. Hence, 

according to Weierstrass’ factor theorem, the first "bracketed term is 
equal to the following expression: 


X+iY 

2rt 


\e 


4- (w - h) + iK' 


,hi(w) 


where hj^(w) is an arh'itraxy entire function. When the other bracketed 
terms are treated in a similar manner 


X+iY 


bW(w) , J 


0 


^(w - h) + iK^ i 


0 


|^^(w - b') + iK-] 


X'+iY' 
2n 




+ b + 2a) +iK* 


h7i(w) 


X-iY 

2rt 


0 


K 


+ b' + 2a) +iK* 


h4.(w) 


X'-iY' 
2rt 


(B20) 


From equation (B20) it follows that 


^ |i°8e 0 


X'+iY' 

2« 

X-iY 
2n 

X'-iY' 

2n 


— (w - b) + iK' 
cD]L 


+ + 


loge©|^^ (w - b') + iK'J + hgCw^ + 
jlogg 0|^^ (w + b + 2a) + iK'j + h 3 (w)| + 
(w + b' + 2a) + iK'j +' h 4 (w)| 


log^ 0 


(B21) 


It is now desirable to find the relation among the four arbitrary func- 
tions hi(w), hgCw), h 3 (w), h 4 (w), and the constant A^. Differ- 

entiation of equation (B21) and use of equations (B17) and (B18) give 
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V(w) = 


JL f 

a>2^ \ 2jt y 


0 




,-f(w-b) +iK' 0' 

-I ( X'+iY^ 

a>2^ V 2jt 


— (w-Ta) +1K’ 


— (w-bO +iK' 

l“i 


t 


0 


^(w-bO +1K' 
o>]_> ■' 


K 

Oi 


i fX'-iY'] 

iV 2n J 


0. 


0 


— (w + b' + 2a) +iK» 


— (w + b! + 2a) +iK' 


0' 

K yX-iY \ _ 

■o>lV 2n J Q 


— (w + b + 2a) +iK? 

L“1 _ 


— (w + b + 2a) +iK! 


W ^ ^ ^ W,(w) 


where Wj^(w) = ^ hjj^(w)j k = 1, . . .,4. 

From reference 4 (p. 480), 

H! fw) in 0! (w + iK' ) 

H(w) “ 2K • 0(w + iK' ) 

so that on comparing equations (B22) and (l3) 

Ai = ^ [h’i(w) (X+iY) ,+ WgWCXt+iY! ) t 

h' 3 (w) (X-iY) + h^(w) (X! -iY! )] 


T + 


(B22) 


(B23) 


NACA ra 2393 


33 


By using equation (B23), equation (B21) can "be rewritten as 
W(w) = W(e») = ^ logg efi (w ■- t) + iKlj + 

W loge e[|. (v+b+2a)+iK'] + log^ (w-V )+iK'] + 
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APPENDIX C 

FORMULATION OF FUNCTION Q(z) 

In order to make certain that the conditions on g(z) (equa- 
tions (32) and (33)) are satisfied, a function Q(z) is introduced 
as follows: 

Let 


g(z) = Q(z) h(z) . (Cl) 

where h(z) is analytic in the annular region and Q(z) has zeros 

at Re^^ and 1. Furthermore, in order that the real and the imaginary- 
parts of h(z) may he separated later, Q(z) must be real on the 
boundary of the annular region. 


Let 


z = e 


and 

Consider the product 

■r 


Q(w) = Q[z(w)] 


^ K \ 
(HI- -)h 


03 , 


(w - r - iy) 


) H 

— (w + 2a) 

CD-, ^ 

H - 
c 

J 

1 

1. 


— (w + 2a + r - i-)f)| 


where (reference 8, ch. XXl) 


H(w) =i^^(w'«^2"^ I T ) 

and 

0(w) = I t) 

From reference 4 (p. 487 ) y 

‘^^(y + z)-di(y - = '^i^(y)'®4^(z) - "34^(y)'&i^(z) 


(C2) 


(C3) 


(C4) 


(C5) 
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so that if 




z = 


-a:r 

2C0, 


Then 


or 


^] (C6) 

ik H “fe " '"H ■" ®1?!] ■ 

H^[^]0^[i(..a,] 


(C7) 


Similarly, 


H 


— (w - r - ir) 
L“l 


H 


(w + 2a + r - ir) 0^ (O) 


CD- 


— (w + a - 2ir) 


] ^ ■■>]- 

+ '>] (w + s - 2ir)] 


(C8) 


so that for z = Ae^® or w = -a + i0 equations (C7) and (C8) are 
• both real . 
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Therefore^ if 


Q(w) = J 

rn 

1 

(V)l 

H 

— (w - r - ir) 

(D-, 

H - 
c 

1 

L 

L 1 J 


^ 1 ■ J 

- 


(w + 2a) 


H 


J L“l 


K 


(w 



q[w(z)] has zero value at z = 1 
values on the circle | z|= A. 

For computation along | ^ | = 

Q(w) = 


H 


K y N „ K / 

— (w) H — (w - 
^1 ^1 


r - ir) 


and z = and takes on real 


C; write Q(w) in the form 


H 


Kl 


(w - 2c) 


H 


K 


(w - 2c + r 


ir)j 

(CIO) 


This equation can be expressed in a form similar to that of equa- 
tion (C9) . 
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Singular ity positions -- First reflection in outer circle 
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(a) Path of integration. 




(b) Profile corresponding to figure 4(a). 

Figure 4. - Illustration of possible lack of closure. 
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